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Abstract 

Wave disturbances of a stratified gas are studied. The description 
is built on a basis of the Bhatnagar - Gross - Krook (BGK) kinetic 
equation which is reduced down the level of fluid mechanics. The 
double momenta set is introduced inside a scheme of iterations of the 
equations operators, dividing the velocity space along and opposite 
gravity field direction. At both half-spaces the local equilibrium is 
supposed. As the result, the momenta system is derived. It reproduce 
Navier-Stokes and Barnett equations at the first and second order in 
high collision frequencies. The homogeneous background limit gives 
the known results obtained by direct kinetics applications by Loyalka 
and Cheng as the recent higher momentum fluid mechanics results of 
Chen, Rao and Spiegel. The ground state declines from exponential 
at the Knudsen regime. The WKB solutions for ultrasound in expo- 
nentially stratified medium are constructed in explicit form, evaluated 
and plotted. 
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1 Introduction. 



Devoted to D.A. Vereshchagin memory 

There are gas dynamics problems at which it is necessary to use a basis 
going out of traditional Navier - Stokes hydrodynamics. It is connected to 
a break of the condition: Kn = l/L « 1, where Kn - Knudsen number, 
I - free particle path, and L - characteristic scale of non-homogeneity of a 
problem. Perhaps, the first work, in which a wave disturbance in a gas was 
investigated from the point of view of more general kinetic approach was 
the work of Wang Chang and Uhlenbeck [I]. The authors have offered a 
method of a dispersion relation construction in a homogeneous gas directly 
from Boltzmann equation. 

The further theoretical and experimental researches [3] - on sound 
propagation in a homogeneous gas have shown, that at Knudsen numbers of 
the unit order the waves behavior considerably differs from ones predicted 
on a basis of Navier - Stokes equations. These researches revealed two es- 
sential features: first, the perturbations keep wave properties at more large 
values of Kn, than it could be assumed on a basis of the classical hydro- 
dynamic description. Secondly, at Kn > 1 such concepts as a wave vector 
and frequency of a wave become ill-determined. May be the most adequate 
results that reproduce experiments [3] almost in all the range were obtained 
in |T2]. It is more difficult to explore the case, when the Knudsen number is 
non-uniform in space or in time and passes the Knudsen regime area. The 
statement and the solution of such problems should definitely be based on a 
kinetic equations or their advance model analogues. 

Quite recently interest to the problems has grown again in connection 
with general fluid mechanics development [THl H3 [TBI EH 1201 ■ ^ was pushed 
by more deep understanding of perturbation theory (so-called nonsingular 
perturbations), see, e.g. [2~T| . 

In his paper we consider the gas medium, stratified exponentially in grav- 
ity field, directed along z axis. In means that the Knudsen number also de- 
pends on z: Kn{z). We continue to develop the method [TSJ that goes up to 
the pioneering paper of Lees [221- The construction of analytical solutions of 
the model kinetic equation Bhatnagar - Gross - Krook (BGK) is extracted 
via separate representation of the distribution function as the local equilib- 
rium one but with different momenta sets at positive and negative velocity 
component v z subspaces. 
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Thus, the set of parameters determining a state of the gas increases twice. 
Such number of parameters of the distribution function (JBJ) results in that 
the distribution deviates from a local-equilibrium and accordingly widen hy- 
drodynamics. In the range of small Knudsen numbers I << L we have 
M+ = M~ and distribution function © passes to local equilibrium one, giv- 
ing a solution of the Navier-Stokes hydrodynamical regime. For big Knudsen 
numbers the formula © gives a solution of so-called collisionless problems. 
Similar ideas have resulted successfully in a series of problems. For exam- 
ple, in papers [22] - |2Z] a method of discontinuous distribution functions 
was used for the description of a flat and cylindrical (neutral and plasma) 
flows (221 " • For a flat problem the surface of break in space of speeds 
was determined by the same natural condition V z = 0, and in a cylindrical 
case V r — 0, where V z and V r , vertical and radial component of speed of 
particles respectively. The problem of a disturbance launched by a pulse 
movement of plane j2Z] was solved similarly. In a problem of a shock wave 
structures [23 [2H1 l2Hj the solution was represented as a combination of two 
locally equilibrium functions, one of which determines the function before 
front of a wave, and another - the tail. In a problem of condensation and 
evaporation of drops of any size [301 EI] a break surface was determined by 
so-called "cone of influence", thus all particles were divided to two types: 
flying " from a drop" and flying " not from a drop" . 

At the first two sections we derive the basic equations using the iterations 
in the evolution operator along the idea of the nonsingular perturbation 
method. Next (Sec. 4) we analyze the transition to a limiting case of a gas 
disturbances at large collision frequencies up to the Barnett case. Next (Sec. 
5) we check the free molecular flow limit, demonstrating the declinations from 
the exponential behavior of the "atmosphere" gas density [321 ?]• At the final 
section we construct solutions of the main momenta system by the method 
VKB, considering the wave scale less than the scale of the inhomogeneity. 

2 Linearized BGK equation 

The kinetic equation with the model collision integral in BGK form looks 
like: 

df -df df 
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here / - distribution function of a gas, t - time, v - velocity of a particle of 
the gas, r - its coordinate vector, 

n I (v-Uf 

fl = i ex P 



J 3 / 2 V 3 \ V 2 

- local-equilibrium distribution function, H - a scale of inhomogeneity (in 



atmospheric models H = kT/mg), v T = y2kT/m - average thermal velocity 
of movement of particles of gas, v = exp(— z/H) - effective frequency of 
collisions between particles of gas at height z. It is supposed, that density 
of gas is denoted as n, its average speed U = (u x , u y , u z ) and temperature T 
are functions of time and coordinates. 

Considering small gas perturbations, we shall approximate a distribution 
function / as: 

f(t,r,v) = f (z,v)(l + cp(t,r,v)) , \ip\ « 1 . 

Here /o - equilibrium Maxwell - Boltzmann distribution function , and ip - 
dimensionless small-amplitude perturbation to equilibrium distribution func- 
tion. For thermodynamic parameters of a gas we shall write: 

n = n (l + n'), U = U'v T , T = T (1 + T'), 

where n 1 , U' and T' - dimensionless small additives to equilibrium density, 
average speed and temperature gas respectively. Expanding the local equi- 
librium function f) up to the first order by small amplitudes and taking into 
account, that equilibrium function f identically satisfies the kinetic equation 
((H) we obtain linearized equation BGK. 

We shall consider one-dimensional case of the equation: 

d i^-^-^{±M M) - V }. ( 2 ) 

Here Xn{v) - eigen functions of the linearized collisions operator: 

2 v z (v 2 5\ 



Xi = 1 , Xi 

X2 : 



V$v T \v 2 2 



[2 (v 2 3\ [Q v z f 2 5 2 

X3 = V3U"2) ' x& = VRT 
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M n - the moments of distribution function determined through scalar prod- 
ucts: 

1 f 

M n (t, Z) =< Xn, ¥ >= nl /2 v 3 J d V eX P (- v2 / V i) ■ Xn{v)¥{t, Z, v) . (4) 

In linear approach the moments M n are linked to the thermodynamic vari- 
ables as: 

M 1= ^, M 2 = V2^, M 3 = M^- (5) 
n v T V 3 T 

3 A method of piecewise continuous distri- 
bution function in linear approximation. 

Following the idea of a method of piecewise continuous distribution functions 
let's search the solution if of the equations © as a combination of two lo- 
cally equilibrium distribution functions, each of which gives the contribution 
to the corresponding area of the velocity subspace as follows: 



n v z . f2(v 2 3\ . 1 ' 

V - = Mf + V2— -M 2 - + W- — -- M 3 ~ , v z <0 



The parameters of locally equilibrium distributions functions are linked 
to the correspondent density, average speed and temperature and three higher 
moments. This idea of the method of two-fold distribution functions of © 
further is realized as follows. Multiplying the BGK equation (J2J) by the 
eigen functions (j3J) with the account of (jlj), calculating integrals we obtain a 
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system of the differential equations for the moments M n : 

T M 2 = , 

2 



2H 



dt 2 dz 
dM 2 d 

dt T dz 
dM 3 v T d . , , , , . v T , , 

-y^(M 2 + M 4 )--M 4 



Ml + M 3 - -M 5 ) - ^ (m 3 - -M 5 ) = , 



at 

<9M 4 

~oT 

dM 5 







+ v T — ~M 3 - -M 5 - -f-M 5 = -v(z)M A , 



9 /5 
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6 
2 



2# 



+ w T — ( -M 2 - -M 5 + -M 6 ) + ( -M 4 - -M, 

(72 



# V5 



-z/(z)M 5 , 



<9M 6 w T <9M 5 , v T 



dt 



(7) 

The moments M n of the distribution function are connected to parameters 
of the two-fold distribution functions (0) by the relations: 



M 2 = v/2 (Mf + ^M+ + 1m 3 ~) 



M 4 
M 5 



Here the following notations are introduced: 
1 



: (M+ + Mi 



Mt = (M+ + M- 
M 3 + = (M 3 + + M 3 - 



-Mf + -M 3 " , 
/6 



£ 2 > 



- ( -M7 + -Mr 
5 V3 1 2 3 



Mr 


= 


-Mf) 


M 2 " 


V07T 


-M 2 ") 


M 3 " 


-M 3 -) 



(8) 
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The feedback between parameters and the moments M n is determined 
from (jSJ) and looks like: 

M+ = M 1 + V3M 5 , Mf = f-M 4 + -V6M e 

10 V o 

2 12 / 10 

M+ = v/2M 2 + -V5M 4 + — W— M 6 , Mf = -v^Mg , 
_ 5 5 V 3 

Mi = (^M 4 + ^M 5 ) , M 3 = ^VE(m a + ^VEM ( 



4 A limiting case of large collision frequen- 
cies. 

Within limits of small Knudsen numbers, from the equations ((Zj) it follows, 
that M 4 , M 5 , M 6 << Mi, M 2 , M 3 . Then in a limit z/ -»• oo (Kn < 1, hydro- 
dynamical limit) we have M 4 , M 5 , M 6 — ► and the system (|7J) tends to the 
linearized Euler's system: 

M lt + ^-M 2z - ^-M 2 = , 

M 2t + v T (Mi + M 3 ) z - ^M 3 = , 

M 3t + yf, = , 

where the bottom indices t, z denote partial derivatives by the time and the 
coordinate z. 

In the next order of the perturbation theory in the small parameter 
(fTo) -1 << 1 (here r = H/v T - characteristic time) we find connections: 

M 4 = -^M 3z , M 5 = ^M 2z , M 6 = , 

substituting which into the system of first three equations (J7J) one arrives at 
the linearized system of Navier — Stokes equations: 

M u + ^M 2z - ^-M 2 = , 

M 2t + v T (Mi + M 3 ) z - ^-M 3 - \ V -^M 2zz = , 

Mat + ^M 2z - \^M 3zz = . 
3 6 v 
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In the higher orders of the theory from the system (J7J) s the linearized 
Barnett's equations follow. For example, in the third order of the small 
parameter (utq)^ 1 << 1 we have: 

M lt + ^-M 2z - |^M 2 = , 

M 2t + v T (Mi + M 3 ) 2 - ^M 3 - \ V ^M 2zz - V \ (m 1z + iMi - ^=M 3 ) = , 

M 3t + ^M 2z - l^M 3zz - % (M 2z - ±=M 2 ) = . 
3 6 v v l V H J zz 



5 A free molecular flow limit. 

Next we would study the opposite limiting case - so-called Knudsen regime 
[y << 1/tq). For simplicity we restrict ourselves now by a consideration of 
stationary solutions of the systems (J7J). 

At v = ( collisionless gas) the general solution of the system (J7J) is 
expressed in elementary functions 

Mi = C 1 + C 3 -i-3C 5 exp(-^) > M 2 = C 2 exp(^ 
M 3 = C 3 -ic 5 exp(-z2iJ) , M 4 = (c 4 - C 2 -|) exp (-|) , 

M 5 = C 5 exp(-^) , M 6 = 1 -{2C, - 7C 2 - 5C 6 ) exp 

Generally, when z/ 7^ 0, the solutions for M 2 and M4 remain the same, while 
for M 5 it is possible to derive the linear non-uniform equation with variable 
coefficients 



1 5 i/q f 2z\ 7i/ 



"2- 



Leaving borders of the system motionless, let's assume that average speed of 
the gas M 2 = 0. Then, if C 2 = 0, we go to the linear homogeneous equation 
for the moment M5. In dimensionless variables z — > z/H it is: 

M 5zz + ^M 5z - Ke~ 2z M 5 = , (9) 
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where 

K=-[^-\ =-(Kn 
3 V v T J 3 V 



5 (v Q H\ 2 5 /tr _ N _ 2 



Let's change the coordinate variable as exp(— 2z) = t, then the equa- 
tion (jHJ) becomes: 

tM 5tt + ^M 5t -kM 5 = 0, (10) 

5 v 2 H 2 



where k , . 

12 i> T ^ 

The equation (fTUj) represents the linear differential equation ( degener- 
ate hypergeometric one) of the second order with variable coefficients. The 
solution may be found as the generalized power series 

oo oo 
n=0 n=0 

Returning to the former dimensionless variable z = — one arrives at: 

H 

oo — oo 

M 5 (z) = ]T a® ex P (-2nz) + exp(--) £ bf exp(-2^) . (11) 

n=0 ^ n=0 

For the coefficients af£' and we have recurrent formulas: 

(5) Ak (s) (5) 4k (B) 

n+1 (n + l)(4n + 3) n ' n+1 (n + l)(4n + 5) n ' 

Substituting M 5 into the equations for other moments yields: 

M 2 (z) = 0, 

M 4 (s) = M 40 exp(z) , 



OO 



M 3 (s) = M 30 - l—zM 40 + \4 ] z + 4«o 5) + E fl i 3) ex P (~2nz) 
5 u T 5 15 ^ 



1,(5) / ^ / ^ X 



- -^ 5) exp (--J + exp (--J ^ &i 3) exp (-2n*) , 

(12) 

where 

(3) _ 1 2n - 3 (5) 3 _ 1 An - 5 (5) 

30 n n ' n 15 4n + l n 



Next the expressions for M 6 and M^. 



5 v H 00 

M 6 (z) = M 60 exp(z)+--?—eMz)Y, b n ) exp (-^(4n + 5)/2) + 

J V T n =0 

+ J^ eX p(^)^ a ( 6 )exp(-2^ + l)) , 

Mriz) = M 10 + ^M 40 (2z-^)+g + ^ 2 -^)a + (13) 
+ 36 exp -- + M m {z - 1) + CL ( n ] exp (-2nz) + 

^ ^' n=l 

+ exp -- ^fo«exp(-2nz) . 

V Z/ n=l 

The coefficients are: 

(1) 1 (2n + l)(6n + l) (5) (1) 1 (12n + 5)(4n + 3) (5) 

20 n 2 n ' n 5 (4n+l) 2 n ' 

„( 5 ) 0/,(5) 



2(n + 1) ' n 4n + 5 ' 

From the recurrent formulas for a£f) and fe^ 5 -* it is possible to derive their 
expressions via constants and 6q 5 . 



(5) = J5) (4fc)» 

ra!-3-7-ll-...(4n + 3) 



6l 5) = &< 5) 



7? 

(Ak) n 



n! • 5 • 9 • 13 • . . . (An + 5) 
n factors 

The series determining the solution (fP2)) and (|13p. converge at any 2 7^ 0. 
Behavior of the moments M n as functions of z are shown in figures Fig 1. 
Some discussion of this stationary case that could be considered as the ground 
state of Knudsen atmosphere theory is published at [32], where a verification 
of the theoretical results are made via molecular dynamics simulations. The 
deviations of exponential behavior of such atmosphere were discussed also in 
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6 Construction of solutions of the momenta 
system by WKB method. 

In this section we apply the method WKB to the system (|7jl. We shall 

assume, that on the bottom boundary at z = a wave with characteristic 

frequency u is generated. Next we choose the frequency u to be large 

3u H 

enough, to put characteristic parameter £ = ^> 1. We shall search for 

v T 

the solution in the form: 

M n = ip n exp(iuj t) + c.c. , (14) 

where, for example, ipi, corresponding to the moment Mi, is given by the 
expansion: 

6 oo i 

^i = EE T-iw^S? «P(«V*(*)) , (15) 

fc=lm=l V 1 ^) 

here <pk{z) - the phase functions corresponding to different roots of dispersion 
relation. For other moments M n , n = 2, . . . ,6 corresponding functions ip n 
are given by similar to (fTSJl expansion. The appropriate coefficients of the 
series we shall designate by corresponding Bffl C$ Dffl Effl Fffl . Substitut- 
ing the series (|T5|) at the system (JZJ)one arrives at algebraic equations for the 
coefficients of (fTK|) in each order. The condition of solutions existence results 
in the mentioned dispersion relation: 



(16) 



6r/ 3 + (5m 2 + 20m - 21)r] 2 + (5m 3 - 24m 2 - 33m + 10)77- 
-3m 3 + 9u 2 + 9iu - 3 = . 

Here for convenience the following designations are entered: 

For the coefficients . . . the algebraic relations are obtained: 

B (k) = ^ fX A (k) ^ c (k) = 1 4m - 1 +Vk A (k) ^ D (k) = ± 1 1 + 577 ft 



3 2iu + 3r]k 2 2iu + 3r]k 

n(k) 1 3m - 5iur] k + 5i]k - 5r/ fc 2 (fc) (fc) i^/2rjk~ (fc) 

= /i^ , Pi = =P Ej-y 

rjk 2iu + 3r/k 3u 



The dispersion relation (jTBjl represents the cubic equation with variable 
coefficients, therefore the exact analytical solution by formula Cardano looks 
very bulky and inconvenient for analysis. We study the behavior of solutions 
at v — ► (free molecular regime) and v — > oo (a hydrodynamical regime). 

At v — the dispersive relation becomes: 

6t] 3 - 21r] 2 + 10r7 — 3 = 0. 

The roots are: 

Vl = 3 , T] 2 « 0.25 + 0.32i , 7? 3 « 0.25 - 0.32i . 

Specifying roots (|16|) by the theory of perturbations up to m 3 for the three 
solutions branches it is obtained: 

m = 3 - 0.56m 2 -i(1.96u- 0.17m 3 )... , 

772 = (0.25 + 0.32i) + (0.07 - 0.69i)u - (0.14 - 0.52z)u 2 + (0. 19 - 0.08i)u 3 . . . , 
r] 3 = (0.25 - 0.32«) - (0.07 + 0.69z)u - (0.14 + 0. 52i)u 2 + (0.19 + 0.08i)u 3 . . . . 

Correspondingly for ip zk = -^-y/rjk we have: 
<p gl = 1. 15- 0.47m 2 -z(0.37m- 0.17m 3 )... , 

ip z2 = (0.38 + 0.19i) - (0.12 + 0.34i)w+ (0.12 + 0.02i)u 2 + (0.06 + 0.07i)w 3 . . . , 
<^ 23 = (0.38 - 0.19i) + (0.12- 0.34i)w + (0.12- 0.02i)w 2 - (0.06 - 0.07i)w 3 . . . . 

Similarly in a limit v — > oo (a hydrodynamical limit) for solutions of the 
equation (fTtjj) 77^ it is derived: 

77! = 0.6 - 1.87m- 2 - i(0.72M- 1 - 4.27m" 3 ) , 
772 = -%{u + 0.80m- 1 ) + 1.40 + 0.18m~ 2 , 
m = -0.83n 2 + 1.50 - i(2.33u - 1.52m- 1 ) . 

The first root relates to the acoustic branch. Accordingly, for the ip z k we 
have: 

(p gl » 0.52 - 0.71m- 2 - ^O.SIm- 1 - 1.41u- 3 ) . . . , 

Vz2 « - i)(0.47 + 0.30m- 2 ) + v/m(1 + i)(0.33M- 1 + 0.07m- 3 ) . . . , 

ip z3 « -0.85 + 0.62m- 2 + i(0.6lM + 0.05m- 2 ) .... 
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The solution of the equation ([16)1 at any u is evaluated numerically. Be- 
havior of real and imaginary parts r) n and (p z k as functions of u, and their 
behavior on a complex plane are presented at the figures 1-4. 




Fig. 1. Behavior of the real part of roots of a dispersion relation: 

a - 7)i, b - 772, c - 1)3 
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Fig. 2. Behavior of the imaginary part of the roots of the dispersive 

relation: 
a - 771, b - 772, c - r? 3 
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Fig. 4. Behavior of the imaginary part of the phase functions 0. 

a - l2 , b - (p 2z , c - (p3z 



As an illustration let us consider a problem of generation and propagation 
of a gas disturbance, by a plane oscillating with a given frequency uj . Let's 
assume, that at height z — all moments = 0, except for M 2 — Uq. In 
evaluations we we choose the dimensionless frequency equal to u(0) = 0.1, 
that fits heights z 300 in the Earth atmosphere. Behavior of the solutions 
for M k constructed by the VKB technique are shown on fig. 7. 
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Fig. 5. Dependence of the moments on height, u(0) = 0.1: 

7 Conclusion 

In this paper we propose a one-dimensional theory of linear disturbances 
in a gas, stratified in gravity field, hence propagating through regions with 
crucially different Kn numbers. The regime of the propagation dramatically 
changes from a typically hydrodynamic to the free-molecular one. We also 
studied three-dimensional case, to be published elsewhere. The theory is 
based on BGK or Gross- Jakson kinetic equation, which solution is built by 
means of locally equilibrium distribution function with different local pa- 
rameters for molecules moving "up" and "down". Equations for six moments 
yields in the closed fluid mechanics system. For the interesting generaliza- 
tions of the foundation of such theory see the recent review of Alexeev |34j . 
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